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Abstract 

We propose an analytical formula for the anomalous scaling exponents of inertial range structure 
functions in incompressible fluid turbulence. The formula is a Knizhnik-Polyakov-Zamolodchikov 
(KPZ)-type relation and is valid in any number of space dimensions. It incorporates intermittency 
in a novel way by dressing the Kolmogorov linear scaling via a coupling to a lognormal random 
geometry. The formula has one real parameter 7 that depends on the number of space dimensions. 
The scaling exponents satisfy the convexity inequality, and the supersonic bound constraint. They 
agree with the experimental and numerical data in two and three space dimensions, and with 
numerical data in four space dimensions. Intermittency increases with 7 , and in the infinite 7 limit 
the scaling exponents approach the value one, as in Burgers turbulence. At large n the nth order 
exponent scales as y/n. We discuss the relation between fluid flows and black hole geometry that 
inspired our proposal. 
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The remarkable phenomenon of fluid turbulence is one of the major unsolved problems of 
physics [l|. Most fluid motions in nature at all scales are turbulent. Aircraft motions, river 
flows, atmospheric phenomena, astrophysical flows and even blood flows are some examples 
of set-ups where turbulent flows occur. Despite centuries of research, we still lack an ana¬ 
lytical description and understanding of fluid flows in the non-linear regime. Insights into 
turbulence hold a key to understanding the principles and dynamics of non-linear systems 
with a large number of strongly interacting degrees of freedom far from equilibrium. In 
addition to being a major challenge to basic science, understanding turbulence is likely to 
have an important impact on diverse practical problems ranging from environmental issues 
such as pollution and concentration of chemicals to cardiovascular physiology. 

In this paper we will mainly consider incompressible fluid flows in d > 2 space dimensions. 
They are the relevant flows when the velocities are much smaller than the speed of sound. 
The incompressible Navier-Stokes (NS) equations provide a mathematical formulation of the 
fluid flow evolution. They read 

dtv^ + v^djV^ =—dip + udjjv'^, = 0, i = l,...,d, (1) 

where n* is the fluid velocity and p is the fluid pressure. An important dimensionless param¬ 
eter in the study of fluid flows is the Reynolds number TZg = where / is a characteristic 
length scale, v is the velocity difference at that scale, and u is the kinematic viscosity. The 
Reynolds number quantihes the relative strength of the non-linear interaction compared to 
the viscous term. When the Reynolds number is of order a thousand or more, one observes 
numerically and experimentally a turbulent structure of the flow. This phenomenological 
observation is general, and fluid details are of no importance. The turbulent velocity held 
exhibits highly complex spatial and temporal structures and appears to be a random pro¬ 
cess. Thus, even though the NS equations are deterministic (in the absence of a random 
force), a single realization of a solution to the NS equations is unpredictable. 

Instead of studying individual solutions to the NS equations, one is led to consider the 
statistics of the solutions. The statistics can be dehned in various ways. One can use an 
ensemble average by averaging over initial conditions. Turbulence that is reached in this way 
is a decaying one. Alternatively, one can introduce a random force. This allows reaching a 
sustained steady state turbulence with an energy source and a viscous sink. The statistical 
properties of turbulent hows are remarkable. Numerical and experimental data show that 
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the statistical average properties exhibit a universal structure shared by all turbulent flows, 
independently of the details of the flow excitations. One defines the inertial range to be 
the range of distance scales Ly r ^ Lp, where the scales Ly and Lp are determined by 
the viscosity and forcing, respectively. Turbulence at the inertial range of scales reaches a 
steady state that exhibits statistical homogeneity and isotropy. 

One defines the longitudinal velocity difference between points separated by a fixed dis¬ 
tance r = |r| 

5v{r) = (n(r, t) — n(0, t)) . (2) 

r 

The structure functions Snir) = {[5v{r)Y) exhibit in the inertial range a scaling Sn{r) ~ 

The exponents are universal, and depend only on the number of space dimensions. 
In 1941 Kolmogorov argued that in three space dimensions the incompressible non- 
relativistic fluid dynamics in the inertial range follows a cascade breaking of large eddies 
to smaller eddies, called a direct cascade, where energy flux is being transferred from large 
eddies to small eddies without dissipation. He further assumed scale invariant statistics, 
that is 


P{5v{r))5v{r) = F 


6v{r) 


( 3 ) 


where P{6v{r)) is the probability density function, and h is a real parameter. Treating the 
mean viscous energy dissipation rate e as a constant in the limit of infinite Reynolds number, 
he deduced a linear scaling of the exponents = n/3. 

All direct cascades are known numerically and experimentally to break scale invariance 
and do not simply follow Kolomogorov scaling. Note, that in two space dimensions the energy 
cascade is inverse, that is the energy flux is instead transferred to large scales. Kolmogorov’s 
assumption that the random velocity held is self-similar is incorrect in direct cascades, but it 
seems to hold in the inverse cascade. The self-similarity assumption misses the intermittency 
of the turbulent flows. Thus, in order to calculate the scaling exponents one has to quantify 
the inertial range intermittency effects. The calculation of the anomalous exponents and 
their deviation from the Kolmogorov scaling is a major open problem. 


Since 1941 many multifractal models of turbulence have been proposed [l|. These express 
multifractality directly in terms of fluctuations of the velocity increments or of the energy 
dissipation. For example, Kolmogorov and Obukhov js, 4| proposed to replace the constant 
global average e with local averages over a volume of dimension r. One then considers as 
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a lognormally distributed random variable with variance ~ —ln{r). Later, Mandelbrot 
argued that one should think about the energy cascade as a random multiplicative process. 
In this case a random measure can be formalised mathematically as a limit of a “Gaussian 
multiplicative cascade” (see joD-The lognormal model assumes “refined self-similarity” 

6v{r) 


P{6v{r))6v{r) = F 


(4) 


.((G)r)V3 

which leads to a formula for producing physical inconsistent supersonic velocities at 
large n and a violation of the convexity inequality [l|. 

Our proposal in this paper is different. We propose that in intermittent turbulence, Kol¬ 
mogorov linear scaling itself is evaluated with respect to a lognormal random measure. This 
is a “gravitational” dressing of Kolmogorov scaling which is inspired by the relation of fluid 
dynamics and black hole horizon dynamics in one higher space dimension [^. We propose 
the dressing of Kolmogorov scaling is via a KPZ (Knizhnik-Polyakov-Zamolodchikov)-type 
relation js]. This gives an analytical formula for the scaling exponents of incompressible 
fluid turbulence in any number of space dimensions d> 2. It reads 

n 


en - 3 = 7Gd)en(l-en) , (5) 

where '^{d) is a numerical real parameter that depends on the number of space dimensions 

d. 

A major part of the paper will be devoted to checking our proposed formula ([5]). We will 
hrst verify that the scaling exponents obtained from (|5]) satisfy the convexity inequality 
and the supersonic bound constraint. We will then show that they agree with the experi¬ 
mental and numerical data in two and three space dimensions, and with the numerical data 
in four space dimensions. Intermittency increases with 7, and in the infinite 7 limit the 
scaling exponents approach the value one, as in Burgers turbulence. At large n the nth 
order exponent scales as y/n. 

We will show that the formula does not apply to the Kraichnan model for passive scalar 


advection by a random velocity field 


10l |. This may have been expected since physics 


of passive scalar turbulence is different from the incompressible Navier-Stokes system. On 
the other hand, it is possible that the more general type of random measure can describe 
passive scalar turbulence. 

The paper is organized as follows. In section II we will explain the coupling to a random 
geometry, discuss the proposed formula and its properties, and perform analytical checks and 
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comparison to experimental and numerical data. While we will establish certain properties 
of the function 7 (d), we will not calculate its precise form in the paper. In section III 
we will apply the formula to the passive scalar model. In section IV we will discuss the 
relation between fluid flows and black hole geometry that inspired our proposal. Section IV 
is devoted to a discussion and open problems. 


I. EXACT FORMULA FOR THE SCALING EXPONENTS 


A. Coupling to a Random Geometry 


Coupling to a random geometry means changing the Euclidean measure dx on a to a 

^ wde. .adon. de.d ,(.) K. covadaace 

(p{x)(j){y) — log \x — y\ when \x — y\ is small (but still in the inertial range). Physically, 

the notion of distance r is modihed in the new measure. Consider a set of scaling exponents 
(Hausdorff dimensions in the mathematical setup) ^0 with respect to the Euclidean measure. 
Denote the same set of exponents, but now with respect to the random measure, by Then 
^ and .^0 are related by the KPZ relation 

e-eo = 7'(d)e(l-0- (6) 


Mathematically, t 
(for a review see 


Tis is a known method to obtain a multifractal structure from a fractal one 
^ and references therein). Our proposal is that one can incorporate the 
effect of intermittency at the inertial of range of scales by coupling to a random geometry 
in this way and evaluating the Kolmogorov linear scaling exponents ^0 = f with respect to 
the random measure. 

Physically, it is highly nontrivial that the steady state statistics of turbulence can be 
viewed as such a combination of the scale invariant statistics and intermittency. Note, that 
intermittent features appear at short length scales, and this is when the effects of the random 
held (j) are prominent. We conjecture that the is proportional to local energy hux 

held e{x), 

" ' ( 7 ) 


e(x) = - {diVj + djViY 


in the direct cascade of the turbulent huid. This has some similarities to the Kolmogorov- 
Obukhov lognormal model. In that case, rehned self-similarity implies the following simple 
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dressing of Kolmogorov scaling 




( 8 ) 


Evaluating the expectation of the lognormal energy dissipation, one hnds 


n 2 ^ 


Q 


3 “ ^ 3 3 


(9) 


As noted in this formula fails as it implies is a decreasing function for large enough 
n, which violates basic physical inequalities. The Kolmogorov-Obukhov formula ([9]) consists 
of the leading terms in the expansion at small intermittency 7 of the KPZ formula ([5]). The 
two formulas agree up to order 7 ^. The KPZ formula is different for high intermittency and 
may be viewed as a completion/generalization of the Kolmogorov-Obukhov formula to the 
strong intermittency regime 3 - 

Here we assume instead that the fluctuating dissipation held e{x), acts as a random 
measure. Let us make a few comments on the mathematical structure of this coupling to a 
random geometry. First, note that there are numerical factors that depend on the number 
of space dimensions, between 7 appearing in the random measure and 7 in ([ 6 ]) 6 |. Since 
what is relevant for us is the formula (Ej), we will keep for simplicity the notation where 7 ^ 
appears in ([ 6 ]). 

The KPZ relation was hrst derived by coupling a two-dimensional conformal held theory 
(CFT) to gravity and analyzing the ehect of quantum gravity on the scaling dimensions of 
the CFT 1^. This has been dubbed “gravitational dressing”. The KPZ relation has been 
generalized in various directions. First, it has been extended to an arbitrary number of 
dimensions without reference to a conformal held theory structure Gj. The KPZ formula is 
then viewed as a relation between a set of Hausdorh dimensions measured with respect to 
a Euclidean (Lebesgue) measure and a random measure. SeconT one can consider a more 
general random held than the lognormally distributed one UJ, ll2| . In this case there is the 
generalized relation 

^0 = ^-log, E[W% (10) 


where E is the expectation and W is the random variable associated with the measure (not 
necessarily a lognormal one). We will not use this generalization in this paper, but it may 
^ Note also, that the KPZ formula can be mapped into the lognormal model by exchanging ^ and ^0 
and multiplying by an overall minus sign. 
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be valuable in the study of steady state statistics of other non-linear dynamical systems out 
of equilibrium. 

We will consider the formula ([2]), where 7 takes values in the range [ 0 , cxd). However, 
when 7 > 1, the mathematical construction of the random measure changes. In the two- 
dimensional quantum gravity language, 7 is related to the central charge of the matter 
system c, and the critical value 7 = 1 is the c = 1 barrier. The regime 7 > 1 is a different 
phase of the theory, dubbed a ” dual phase”. There may be a duality relation between two 
phases parametrized by 7 and 7 ' that satisfy 77 ' = 1. This could have an interesting impact 
on the study of turbulence in diverse dimensions. 


B. An Exact Formula 


We propose that the scaling exponents of incompressible fluid turbulence in any number 
of space dimensions d satisfy the KPZ-type relation ([S]). Solving for we get 

((1 + 72)2+ 472(1-1))^+72-1 


en = 


272 


( 11 ) 


where in choosing the branch we required hnite exponents ^n- 7 (d) is a numerical real 
parameter that depends on the number of space dimensions d. It can be determined from 
any moment, for instance, from the energy spectrum. 

There are several immediate properties of the formula flTT]) that we can see. First, using 


u = 3 in ffTTl) one gets the exponent ,^3 = 1 in any dimension, an exact resu 


Kolmogorov which agrees with numerical simulations and experiments. In 13(] this scaling 


t derived by 


was derived without employing the cascade picture, but simply from the fact that the NS 
equations are conservation laws. Second, the scaling exponent .^2 is a monotonically increas¬ 
ing function of 7 , while the exponents n > 3 are monotonically decreasing functions of 
7 . Third, in the limit u —)■ 0 we get that 0, as expected. Fourth, in the limit 7 —?■ 0 
we have ^ that is scale invariant statistics with no intermittency. Fifth, in the limit 
7 —>■ 00 , we have 1, as in Burgers turbulence. The scaling exponents take values in 

the range | < ^2 < 1, and 1 < < f for n > 3. We will propose that the limit 7 —)■ cxo, is 

the limit of inhnite number of space dimensions d. The subleading correction, relevant for 
developing a systematic ^ expansion reads 

1 /n \ 


7 


1 

^4 


( 12 ) 







Sixth, in the limit n —?■ oo for hxed 7 , we have 


1 

n — 

7 





thus growing as y/n. Seventh, at the ’’critical” value 7 = 1 we get ■Cn = (f) 


( 13 ) 


C. Analytical Constraints on the Scaling Exponents 

If there exist two consecutive even numbers 2n and 2n + 2 such that ^ 2 n > Gn+ 2 , then 
the velocity of the flow cannot be bounded. Using flTT]) it is straightforward to show that 
^ 2 n < ^ 2 n +2 for any 7 , thus flTT]) satishes the absence of supersonic velocity requirement. 
The second condition is that of convexity. For any three positive integers < 77.2 < u-s, the 
scaling exponents satisfy the convexity inequality that follows from Holder inequality 

(ns - ni)^ 2 n 2 > (’^a - 772 ) 6 ni + (772 - 77 i)^ 2 n 3 - ( 14 ) 

Using flTT]) it is straightforward to show that the Holder inequality holds. Equality is achieved 
when 7 = 0 , when 7 —)■ cxo and when n^ = rij for some i ^ j and arbitrary 7 . 


D. The Energy Spectrum 

The structure function S' 2 (r) ~ gives the energy spectrum of the fluid. Using flTT]) we 
see that .^2 is a monotonic function of 7 that takes values in the range | < >^2 ^ 1 when 7 
goes from zero to inhnity. In momentum space a deviation from the Kolmogorov spectrum 
for small 7 (small d) reads 

E{k) ~ . (15) 

For large 7 (large d) we have 

E{k) ~ . (16) 


E. Comparison to Experimental and Numerical Data 

The anomalous scaling exponents flTT]) depend on the parameter 7 , which is a function 
of d. We do not know the exact expression of 7 , but it can be calculated knowing one of 
the structure functions, such as the energy spectrum 7 = ( ) • With this knowledge 
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we can then make an infinite number of predictions. In the following we will compare 
the analytical expression flTT]) to the available numerical and experimental data in various 
dimensions. 


1. Two Space Dimensions 

In two space dimensions the energy cascade is an inverse cascade, where the energy flux 
flows to scales larger than the injection scale. In this case, one has the energy spectrum 
agreeing with the Kolmogorov scaling ^2 = f- Using flTTD . this implies that 7 ( 2 ) = 0, and 
that all the other scaling exponents follow the Kolmogorov scaling = f • 


2. Three Space Dimensions 

In three space dimensions we first use the data for the anomalous scaling exponents 
quoted in |l^ from wind tunnel experiments at Reynolds number ~ 10^. This experimental 
data is consistent with numerical data from simulations of the Navier-Stokes equations, see 
e.g. [^. We fit ([5]) to this data using a least squares fit with the function FindFit in 
Mathematica. We see in Figure 1 an excellent agreement, finding that 7 ^ is about 0.161. 

Next, we consider the numerical results for low order structure function exponents and 
non-integer n given in 16|]. The numerical data is consistent with experiment at Reynolds 
number 10^. For this data, the fitted value of 7 ^ is about 0.159 and in Figure 2 we again 
see excellent agreement. Note that if our conjectured relation between the random measure 
and the local energy dissipation field is correct, one can determine 7 ^ independently by 
measuring the scaling exponent of the two point function, (e(a;)e(0)) ~ . This value has 

been found to be ~ 0 . 2 , which in our formula is still consistent with the data. 

3. Four Space Dimensions 


In four space dimensions, numerical simulations of the Navier-Stokes equations were 
performed in [^. The authors found an increase in intermittency, i.e. for n < 3, 

while for n > 3. We took the data for the structure function exponents in 4d 

given in 0 and performed a fit to ([5]). This is shown in Figure 3. Although taken at a 
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FIG. 1. Fit of ([5]) (blue) to experimental data {l^ (Table 2). The dashed line represents Kolmogorov 
scaling. The best fit value of the free parameter 7 ^ is about 0.161. The error on the data is about 
±1 percent. 





FIG. 2. Fit of ([ 5 ]) (blue) to numerical data of low moments m (red). The dashed line represents 
Kolmogorov scaling. The best fit value of the free parameter 7 ^ is about 0.159. 
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FIG. 3. Fit of jS]) to the 4d exponents given in [l7| . The solid line is the 4d fit with 7 ^ about 
0.278. 

relatively low Reynolds number, the results are in agreement with a simple increase in the 
7 ^ parameter in our formula (I5]). The value of 7 ^ in four space dimensions is htted to about 
0.278. Note that their numerical data for same simulation in three space dimensions predicts 
7 ^ about 0.188, which is higher than the experimental data above. This could be related to 
the relatively low Reynolds numbers involved. 


4 . Intermittency and the Large d Limit 

In order to observe intermittency one has to study the short distance statistical properties 
of the fluid flow. There are various measures for intermittency, such as Fn{T) = ^ ^ 

3. Fn{r) are expected to grow as a power-law in the limit r —)■ 0, while staying in the inertial 
range of scales. 

We can analyze the properties of Fn{r) using flTTll . They scale as ~ r", where a is a 
decreasing function of 7 . In the limit 7 —)■ 0 one has a —?■ 0 and no intermittency, while 
as 7 —)■ cxo we get the maximal intermittency a = Numerically, one sees in a 

clear growth of Fn{r),n > 4 in the limit r ^ 0, when as we increase the number of space 
dimensions in the simulation. The data is not accurate enough to observe the growth when 
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n = 3. 

Another exponent that is used to quantify the intermittency is 

A = 2 - & . 


( 17 ) 


Experimentally in three space dimensions it has been measured in the range 0.2 to 0.25. 
Using flTT]) with 7 ^ = 0.161 we get p = 0 . 222 . Expanding around 7 = 0 (d = 2) we have 
/i = 27 ^ + 0 ( 7 ^), while expanding around inhnite 7 we have fi = 1 — ^ 

(also see [^) it was conjectured that in the limit of inhnite d all the exponents 


In 


:7i+ 

m 


approach the same value, one, as in Burgers turbulence [l 8 [. With our formula (|5]) this 
means that 7 goes to inhnity in the limit of inhnite d, and therefore = 1 for any n. This 
suggests the interesting possibility of having a systematic ^ expansion ([ 12 ]). 


II. PASSIVE SCALAR TURBULENCE 


It is natural to ask whether our proposed exact formula for the scaling exponents of 
incompressible huid turbulence is applicable for other systems that exhibit turbulent struc¬ 
ture. In the following we will consider the Kraichnan model for passive scalar advection by 


a random Gaussian held of velocities n*, which is white-in-time {q, 110|. The statistics of the 
velocities is determined by a zero mean v^{t,r) = 0 , and by the covariance 


(n*(t, r)n*(t, T')) = 6{t — t')D^^ . 


(18) 


In the inertial range ~ |rj^, where ( takes values between 0 and 2. 

Examples of passive scalar systems are smoke in the air, salinity in the water and tem¬ 
perature when one can neglect thermal convection. The evolution equation describes a 
passively-advected scalar held T driven by the velocity held 


dtT + v^diT = ndjjT + f , 


(19) 


where k is the molecular dihusivity of T and / is an external force. 

One dehnes the dimensionless Peclet number Ve as the ratio of the scale of huctuations 
of T produced / and the dihusion scale. When Ve ^ ^ there is a scalar turbulence with 
a scalar cascade and constant hux of T^. Similarly to the incompressible huid turbulence, 
one is interested here in the stationary statistics and the scaling properties of the scalar 
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structure functions in the inertial range of scales. Define 6 T{r) = (T{f,t) — T{0,t)) as the 
difference between the values of the scalar held at two points separated by a hxed distance 
r = \ f\. Then, 

Sn{r) = {{ST{r)f^) ~ . (20) 

Here the scale invariant statistics is Gaussian with ^ 2 n = n{2 — (). We can now attempt 
to include the intermittency by the random geometry dressing and the KPZ-type equation 

- n{2 - C) = - ^2n) • (21) 


Solving for ^ 2 n we get 


^2n 


((1 + 7 ^)^+ 47 ^( n ( 2 - 0 - l ))^+ 7^-1 
272 


( 22 ) 


This formula is similar, but not exactly the one proposed by Kraichnan 
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n —)• cx) for hxed 7 , we have ^ 2 n growing as y/n. This is not the expected behaviour, rather 
^ 2 n should approach a constant. This is due to presence of solutions with fronts and shock 
behaviour, which seems to be characteristic of compressible huid systems [^. While the 
KPZ formula does not apply to the passive scalar, it is possible that intermittency in this 
case ultimately can be described by other, more general types of random geometry, e.g. via 
the formula ffTOj) above. 


III. BLACK HOLE HORIZON DYNAMICS 

In the following we will briehy review the relation between huid hows and black hole 
horizon geometry (for a review, see and references therein), that inspired our proposal to 
incorporate the intermittency at the inertial range of scales by a gravitational dressing using 
a random geometry. Consider the Einstein equations with a negative cosmological constant 
A in (d + 2) space-time dimensions 

Eab = Rab — -QabR + RgAB = 0, A, B = 0, ...,d + 1 , (23) 

where qab is the Lorentzian metric, Rab the Ricci curvature and R = q^^Rab the Ricci 
scalar. Black holes are a classical solution of Einstein equations fl23ll . and their hallmark is 
the existence of a horizon H. For example, a black hole solution to the Einstein equations 
has the form 
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( 24 ) 


d 

ds^ = gAsdX^dX^ = —r‘^f{r)dt^ + 2 dtdr + ^ dxidx^ , 

i=l 

where the coordinates are X^ = (t,r,x®). The function /(r) = 1 — (The 
horizon is defined as the surface where /(r) vanishes. It is a (d + l)-diniensional null 
hypersurface, forming a causal boundary preventing any light and particles that cross it 
from returning. Hence, it effectively introduces dissipation. One can associate with the 
black hole horizon a temperature T (appearing in f{r)), and an entropy proportional to its 
cross-sectional area A. In Planck units h = Gn = c = 1, the relation between the area and 
the entropy is Sbh = j- This structure is called black hole thermodynamics. 

Black hole hydrodynamics is a generalization of black hole thermodynamics, similar to 
the generalization of field theory thermodynamics to hydrodynamics. While black hole 
thermodynamics quantifies the thermal equilibrium situation, black hole hydrodynamics 
describes slow derivations from equilibrium. In particular, one can allow the black hole 
temperature to be a slowly varying function T(t, x®) = const.(l+p(t, x®)), and consider black 
hole itself to be moving at velocity x®(t,x) with respect to some rest frame. The perturbed 
solution to the Einstein equations in this setting will yield a slowly evolving curved geometry, 
with the gravity variables providing a geometrical framework for studying the dynamics of 
fluids. This can be made precise in the context of a holographic correspondence, where the 
fluid system lives on a (d -|- 1) dimensional surface of r = const, in the {d + 2) dimensional 
bulk solution. 

The motion of fluids translates to the evolution of the black hole horizon, and the fluid 
variables to its geometrical data. The normal vector n to the horizon H satisfies = 

0. The horizon hypersurface is defined hy r = th = const, and we have 

n®” = 0, n* = 1, n® = X® . (25) 

Thus there is a geometrical representation of the fluid velocity in terms of the normal 
to the black hole horizon hypersurface. As we perturb the black hole and get it out of 
equilibrium, the horizon location changes, and up to an overall constant can be parametrized 
by th 1 +p(t,x). The variable p(t,x) quantifies the deviation from equilibrium and is 
identified as the fluid pressure. The deviation of the horizon area measure from equilibrium 
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is parametrized by the pressure. One gets that the horizon area measure ^ scales like 
dp{t,x^), where d is the number of space dimensions. 

The set of the Einstein equations projected on the horizon 

Eabu^ = 0 , (26) 


describe the evolution of the perturbed horizon geometry, and are equivalent to the incom¬ 


pressible NS equations 


20 


211]. Note that since the Einstein equations are relativistic, 


this amounts to takin g a p articular non-relativistic scaling limit of the fully relativistic fluid- 


gravity equations 


22 


23(1 . For example, one of the equations enforces the vanishing of the 


fractional rate of change of horizon area at lowest order and reduces to the incompressibility 
condition diVi = 0. 

In the gravitational framework, every fluid conhguration that solves the incompressible 
NS equations corresponds to a particular dual bulk solution and corresponding horizon 
hypersurface geometry. To describe forced turbulence holographically, one can introduce a 
background matter term Tab into the Einstein held equations (l23il |24| . This can be used 
to model the large scale random force pumping energy into the huid. The energy injected 
by the pumping ultimately falls into the black hole, where it is dissipated as heat. In the 
turbulent inertial range, the pressure and velocity are stochastic helds, implying that the dual 
horizon geometry and measure itself are also random. Therefore, the statistical properties 
of the random horizon hypersurface (characterized by a sum over surfaces) may encode 
the universal statistical structure of turbulence. At the horizon, the energy dissipation is 
related to the flux of matter across the surface, [25|. Since the horizon 

itself is random, we were inspired to introduce the random measure in order to quantify the 
intermittency effects. 


IV. DISCUSSION 

We proposed an analytical formula for the scaling exponents of inertial range incom¬ 
pressible fluid turbulence in any number of space dimensions d > 2. The idea is that 
intermittency can be taken into account by a novel gravitational dressing of the scale in¬ 
variant Kolmogorov spectrum. Mathematically, the coupling to a random geometry with a 
random measure based on a log correlated held, maps the fractal structure of the scaling 
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exponents to a multifractal one. 

There is one parameter that depends on the number of dimensions that we denoted by 
'y{d). It can be deduced knowing one moment, for instance the energy spectrum. With this 
knowledge one can make inhnite number of predictions. 


Our formula passes the standard analytical consistency checks, such as the convexity 
inequality and the absence of a supersonic mode. Its predictions agree with experimental and 
numerical data in two, three and four space dimensions. The main challenge is to identify the 
origin of the random geometry in a more precise way and determine analytically the function 
'y{d). We expect the holographic framework that inspired our formula - the relation between 
fluid dynamics and black hole horizon geometry - to provide a clean calculational scheme. 
One may also try to calculate 7 using some physical models for the anomalous scaling, such 
as contributions from vortex hlaments 




26| . or statistical conservation laws jl9|. 


While equilibrium statistics is characterized by the Gibbs measure, there is yet no analog 
of this for non-equilibrium steady state statistics. We speculate that there is a general 
principle that allows us to consider the steady state statistics of out of equilibrium systems 
as a gravitationally dressed scale invariant one. If correct, this will shed much light on out 
of equilibrium dynamics. 


Finally, it will be interesting to use the gravitational dressing to study the intermittency 
effects on the anomalous scaling of the transverse structure functions and multipoint cor¬ 
relation functions. Also, our proposed formula is valid for the inertial range of scales, and 
most likely does not incorporate statistical signatures of the dissipation range of scales. It 
is of interest to know whether the latter can be parametrized by a random geometry, since 
after all the Reynolds number is finite in nature. 
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